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> ' Abstract 
00 '. 

QQ ! The problem of the hghtest Higgs boson mass in the next-to-minimal supersymmetric 

] standard model (NMSSM) is investigated. Assuming the validity of the perturbation 

^ I theory up to unification scales and using the recent experimental results for the top quark 

0\ • mass, the restrictions on the NMSSM coupling constants are obtained. These restrictions 



are used to make the predictions for the lightest Higgs mass, which are compared to those 
of the minimal supersymmetric standard model (MSSM). 



1. The aim of this letter is to investigate the hghtest Higgs boson mass problem in 
the next-to- minimal supersymmetric standard model (NMSSM). This model contains an 
additional Higgs singlet, as compared to the minimal supersymmetric standard model 
(MSSM). The cause of consideration of the NMSSM is connected in particular with desire 
to avoid an explicit mass term for the Higgs doublets in a superpotential [|l|. There is 
also an interesting possibility to realize the scenario with spontaneous CP-violation in 
this model 0-0. It is known that in the minimal supersymmetric standard model such 
a scenario [0 contradicts with existing experimental constraints on the MSSM neutral 
Higgs masses p, ^, while in the NMSSM this contradiction can be avoided in different 
ways ^. The experimental constraints on Higgs masses play important role also for 
the case of the absence of spontaneous CP violation. The recent LEP experiments give 
the lower bound on the lightest MSSM Higgs mass rrih > 75GeV |^. Depending on the 
Higgs vacuum expectation values (vev's) ratio tan/3, this bound can be much stronger. 
Thus, one derives [0 m/j > 88GeV for tan/? < 2. This makes the MSSM close to being 
excluded, when the low tan j3 scenario is considered. It is interesting to analyze how the 
situation is looked in the NMSSM. 



2. The next-to-minimal supersymmetric standard model had been widely studied in 
literature lli-|20l. It was noted already in ref. [l^l that the upper bound on the lightest 



NMSSM Higgs particle mass can be different from the one in the MSSM. On the other 
hand, the NMSSM Higgs sector is described (even at tree level) by large enough number 
of unknown parameters, which makes the analysis of Higgs masses very difficult. Today 
only the case, when the supersymmetry breaking parameters of the theory are related by 
universality conditions, has been investigated in details |I7|, ITH]. Requiring for the 



physical minimum of the scalar potential to be the global one and taking into account 
the experimental bounds on the NMSSM particles masses, it was found that the Higgs 
singlet sector is decoupled and - due to smallness of the relevant coupling constants - the 
predictions for the detectable (nonsinglet) Higgs bosons are almost the same as those of 
the MSSM. However, as it was recently pointed out |jl9[, the NMSSM parameters space, 
where the scenario with the universal supersymmetry breaking can be realized, is strongly 
restricted, if it is required for the physical minimum of the scalar potential to be the global 
one. One can avoid such a strong constraints on parameters of theory either assuming 
that the physical minimum of the potential is a local one, or relaxing the universality 
conditions at high energy scales. 

In this letter no universality conditions on the supersymmetry breaking parameters are 
assumed. In this case either only some special regions of parameters space of theory 
are investigated or an upper bound on the mass of the lightest Higgs boson, as a func- 
tion of some parameter(s) of theory, is only found. Here the latter strategy is used: an 
upper bound on the lightest Higgs mass, as a function of tan/3, is obtained. The only 
restrictions on parameters of theory, which are taken into account, are constraints on the 
coupling constants, coming from the assumption of the absence of new physics between 
the supersymmetry breaking scale (Msusy) and unification scales (Mg), as well as from 
the recent experimental results for the top quark mass. The effects, connected with new 
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experimental results for the top quark mass, play important role, when being taken into 
account appropriately. These effects together with two-loop order corrections to the light- 
est Higgs mass lead to the predictions for this particle mass, which are different by about 
(20-25) GeV from those, reported in literature previously. The derived upper bound on 
the lightest Higgs mass is compared to the one of the MSSM. 

3. The NMSSM superpotential is obtained from the MSSM one by making in the latter 
one the following replacement: 

HH,H2 ^ AiV^i^2 + ^iV^ (1) 

(here Hi and H2 are the Higgs doublets and N is the Higgs singlet superfields with scalar 
components Hi, H2, N). Thus, one avoids an explicit mass term in the superpotential 
- such a term was needed, to provide for spontaneous breaking of the SU{2)l x U{1)y 
symmetry. Instead one introduces the additional complex degree of freedom, connected 
with the SU{2)l x f/(l)y Higgs singlet N, as well as two Yukawa-type coupling constants, 
A and k. These couplings enter the Higgs bosons mass matrices in nontrivial way and 
to find the restrictions on the lightest Higgs mass, one has to obtain some constraints 
on A and k and/or correlations between them. Such constraints can be obtained from 
the NMSSM renormalization group equations (RGE's) analysis. Such an analysis in the 
standard model and its extensions has been performed in ref's ||21|]-[Q - there various 
approaches had been developed. 

In this letter it is assumed that all the (gauge and Yukawa-type) couplings are small be- 
tween the electroweak breaking scale ~ lOOGeV (this scale is identified also with Z boson 
mass Mz or with the top quark mass rrit) and unification scales {Mq ~ (10^^ — 10^'^)GeV), 
so that the perturbation theory is applied. This condition leads to restrictions on coupling 
constants of theory. As a result, one obtains some bounds on the lightest Higgs boson 
mass. 

During the renormalization group equations analysis I neglect the supersymmetric parti- 
cles mass thresholds effect, assuming that supersymmetric RGE's are valid from the elec- 
troweak breaking scale up to unification scales. I perform the analysis at one-loop level - 
as the numerical analysis shows, the two-loop corrections do not exceed few percents, i.e. 
are of the same order, as the inaccuracy of approach, connected with the uncertainty of 
choice of the cut-off scale (whether it is chosen ~ lO^^GeV or, say, ~ lO^^GeV). 



4. Let us discuss briefly the problem of the lightest Higgs boson mass in the MSSM 
and its simplest extension with the additional Higgs singlet. In the minimal supersym- 
metric standard model |3l|, |3^ there are five physical Higgs states: two CP-even, one 
CP-odd and one complex charged Higgs states. At tree level these particles masses are 
described by only two unknown parameters: the CP-odd Higgs mass and the Higgs 
doublets vev's ratio tan/3 = V2/V1. As a result, Higgs boson masses and the restrictions 
on them are given by compact expressions. In particular, for the lightest (CP-even) Higgs 
mass it is obtained that rrih < Mz \ cos2/3| < Mz- Notice however that this relation is 
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violated, when radiative corrections to the Higgs potential are taken into account. Usually 
one takes into account only the radiative corrections, connected with the top and bottom 
quarks and squarks loops - this is because of largeness of the top and (for large tan /?) 
bottom Yukawa couplings, as compared to other coupling constants. When being taken 
into account at one- loop level these corrections can increase the lightest Higgs 

mass about (30-60) GeV, depending on parameters space of theory. Recently the two-loop 
order corrections to the lightest Higgs mass had been also calculated [^-[^Ol- 



found that these corrections are important too: they can lower ruh more than lOGeV. 
The above-mentioned radiative corrections to the lightest Higgs mass depend on the fol- 
lowing set of parameters: the top mass rrit, the stop masses m^^, m^^, the sum of Higgs 
doublets vev's squared r]"^ = vf + V2 = {ITAGeVy, the Higgs vev's ratio tan/3, the QCD 
coupling constant = g^/{47T), as well as (via the stops mixing effects) on the stop 
trilinear supersymmetry breaking parameter At and the explicit mass parameter /i. In 
this paper it is assumed that m?^ — m?^ -C m?^ + m?_^ (then stops masses are identified 
usually with the supersymmetry breaking scale). In this approximation one obtains at 



two-loop level m, 40 



rf ml ) 47r^ r/^ 
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+ log + -—: - 327r«3 Xt log + log2 ^^^^ 



mt 

where 



(2) 



_ 2{At-ficotp)^ ( {At- ^icoi(3f \ 
- 772 1 T7T72 l3j 

^^^SUSY \ J-ZiW^.^gy / 

It is easy to see that < 6 - the maximum value of Xt is obtained, when At — fi cot (3 = 
\/QMsusy- In this paper two values of Xt are considered: = - the case of so-called 
no-mixing, and Xt = Q - the case of so-called maximal mixing. 

In the next-to-minimal supersymmetric standard model there are additional one CP- 
even and one CP-odd Higgs degrees of freedom, as compared to the MSSM. At tree 
level the 3x3 symmetric mass matrix for CP-even Higgs fields $1, $2 and Ni [[10| -||13||, 
[pO| is given by formula (Al) of Appendix. Usually one takes $1 = ■^{Re{Hi) — vi) 



$2 = -^{R^{H2) — V2) and Ni = -^{Re{N) — V3), where V3 is the singlet vev. In this 
paper a little bit different representation is used, namely, 

$1 = -^{Re{H^) cos P + Re{H^) sin p-ri) 

$2 = -^{-Re{H^) sin P + Re{H^) cos P) (4) 

Xi = -^{Re{N)-v,) 

Such a choice of the representation of CP-even Higgs fields makes the matrix (Al) more 
transparent for the qualitative analysis. 
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At tree level the upper bound on the lightest (CP-even) Higgs boson mass is the following: 

ml < Ml cos^ 2(3 + X^rj^ sin^ 2/3 (5) 

where the right hand side of eq. (5) is the maximum value of the lowest eigenvalue of 
2x2 upper block (M|. ., i, j = 1, 2) of the matrix (Al) ||Tl], |l^. In other words, the upper 
bound (5) is saturated, when the singlet sector is decoupled^ and the detectable Higgs 
bosons masses are described by the 2x2 MSSM upper block of the matrix (Al). Like 
in the MSSM, there is only one (nonsinglet) Higgs boson with the mass of the order of 
electroweak breaking scale, when the upper bound on nih is saturated]^. 
Even in the region of parameters space, where the singlet sector is decoupled, there is 
a difference between the radiative corrections to the lightest Higgs mass in the MSSM 
and those in the NMSSM. Such a difference arises due to the loops, connected with the 
term (A^ — gl/2)\HiH2\^ of the Higgs potential [0. The contribution of these loops is 



proportional to A"^, X^gl, X'^gf, X'^hf, etc. ((72 and gi are respectively SU{2)l and f/(l)y 
gauge coupling constants, ht is the top quark Yukawa coupling). 

In the next section it is obtained that A is of the order of weak coupling constant g2 
or smaller. This allows one to neglect all aforementioned corrections, except of those, 
proportional to X^h\. These latter corrections arise due to the top quark and squark 
loops. The contribution of these loops is generalized in the case of the NMSSM in the 
following way. 

As it was mentioned above, both in the MSSM and in the NMSSM the upper bound on 
the lightest Higgs mass is saturated, when there is only one light (detectable) Higgs boson 
with the mass of the order of electroweak breaking scale. Then at tree level one may write 
ml < ml^^^ = AXsurf, whereQ Xsm = M| cos^ 2/5/(4r/2) and Xsm = Mlcos^2f3/{Ar]^) + 
X^ sin^ 2/3/4 for the MSSM and the NMSSM respectively. As for the radiative corrections, 
they can be subdivided in following three parts: 

a) the corrections, coming from the one-loop renormalization of Xsm, 

b) one- loop corrections, connected with the stops mixing effect, 

c) next-to-leading order (two-loop order) corrections. 



The corrections a) can be read off from the standard model RGE's 291: one obtains that 



^SM - Xsm (1 - ^/^L. log ^) + TT^/^t.. ^^g ^ 
1 Stt^ mf I loTT^ mf 



^The necessary condition for this to occur is -Mj^^ 0. If even Ms22^-^'s33 ^ -^ln ~ i ^'^^ 
derives ml = Mj^^ — (Ms^J"^ /Mg^_^, so that the upper bound (5) is saturated, when 0. This 

occurs if either A 0, or 113 — > (this case is disfavoured by the experimental bounds on the hghtest 
chargino mass) or in the case of fine-tuning of the relevant NMSSM mass parameters. In other words, 
the upper bound (5) is saturated only in narrow region of the NMSSM parameters space. It is worth to 
stress also that this is the upper bound on the lightest detectable (nonsinglet) Higgs boson. 

^Generally speaking, this rule is violated for tan/3 ^ 1: then the upper bound on rrih can be saturated 
also when more than one Higgs bosons with the masses ~ Mz exist. However these values of tan/3 are 
out of the interest: as it is easy to understand, the difference between the models predictions occurs due 
to the coupling A, which enters the above-mentioned 2x2 block with the factor sin 2/3. 

•^The choice of the coefficient in front of Xsm depends on the choice of the corresponding coefficient 
of the standard model Higgs potential. Here the same notation for Xsm, as in ref. [E9|, is used. 
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where htgj^j = mt/r]. The corrections b) and c) depend only on the top and stop sectors 
of Lagrangians of the models. These sectors (and hence the corrections b) and c)) are 
almost the same in both of models. The difference is only that the mass parameter /i 
is replaced by the product AA^ (or vice versa). Respectively, in equation (3) one has to 
make the replacement /i —>■ Xv^. 

Using the above-mentioned analogy between the models, one obtains the following bound 
on the lightest NMSSM Higgs particle: 



ml < (M| cos^ 2(3 + XY sm' 2/3) f 1 - + Al^l 



^2 ^2 ^2 j ^^2 ^2 



2 

+ 'og ^ + (14 - 32^3) ( X, log ^ + log' Ml^) 

mf loTT^ \1 7]'^ I \ mf J 



(6) 



One can see from comparison of (6) to (2) that the difference between the models pre- 
dictions can occur due to the term, proportional to X^r]"^ sin^ 2/3. This term is important 
for low values of tan/?. Therefore I restrict myself by consideration of tan/3 < 6. To 
estimate the difference between the models predictions, one must find some constraints 
on the coupling constant A. As it was mentioned above, such constraints are derived from 
the renormalization group equations analysis. 

One assumes often (see e.g. |2^) that |A| < 0.87- the infrared (IR) quasi-fixed point, 
derived in the limit ht = 0. Such an assumption is not quite correct: the restrictions on A 



become stronger, when ht ^ 1 [|TT], [T^. As it is shown in the next section, the restrictions 
on A for the experimentally allowed range of the top mass are about 20% different from 
the bound, reported above. 

5. For the low tan/3 scenario the RGB's for SU{3)c, SU{2)l, U{1)y couplings, g^, §2 
and Qi respectively, and Yukawa- type couplings /i^, A and k are the following |3y, ^: 



at 

levr^^ = X{2K' + AX' + ?,hl-\gl-?,gl) (7) 
at 5 

IGtt^^ = ht{Qhl + X^-^^gl-:igl-^^gl) 



, p 2 '^di 3 



where Ci = —33/5, C2 = — 1, C3 = 3 and t = 1/2 log((5^/M|) (as it was mentioned above, I 
neglect the supersymmetric particles mass thresholds effect, considering these RGB's valid 
from the electroweak breaking scale up to unification scales). The behaviour of gauge 
couplings is determined by their experimental values at the electroweak breaking scale 
[0: gi{mt) ~ 0.46, g2{jnt) ~ 0.65, gs^mt) ~ 1.22. Some restrictions on the top quark 
Yukawa coupling is found from the experimental constraints on the top mass. According 
the recent experimental data |42|, mf"'*^ = (173.8 ± 5.2)GeV. One has to transform the 
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pole top mass into the on-shell top mass 



nit = rntirrit) = ht{mt)r] sin P (8) 
to make the predictions for ht. This transformation is done, using the well-known relation 



between the pole and the on-shell top masses [43 



pole 

m{mt) = ^^^^ = (165 ± 5)GeV (9) 
where in eq. (9) only the leading order QCD gluon corrections are taken into account. 



Higher order QCD corrections as well as the stop/gluino loops corrections [^4|, ^ modify 
the on-shell top mass by about 2GeV - they can also cancel each other. Here I neglect 
these corrections: to take into account the stop/gluino loops effects, one has to determine 
both the stops masses and the stops mixing angle. This deviates from the framework of 
this letter. 

As it follows from equations (8), (9), 0.92/ sin/3 < htirrit) < max {ht^^^, 0.98/ sin/5), where 
^Unax is the triviality bound on ht{'mt). This bound, as well as bounds on A and k {Xmax 
and nmax) are found from the analysis of RGE's for these coupling constants. 
The procedure is the following. One evolutes (numerically) A, k, ht and gi, g2, gs according 
their renormalization group equations from the electroweak breaking scale up to the GUT 
scale ~ lO^^GeV. It is required at whole considered energy range for the couplings to be 
small enough to perturbation theory being applied. More precisely, the following condition 
on A, K and ht must be satisfied: 

A'(Q') < 47r, k\Q^) < 47r, /^^(Q^) < 4^ (10) 

The conditions (10) put some constraints on coupling constants at the electroweak break- 
ing scale. Thus, one obtains ht{mt) < 1.15 and \n{mt)\ < 0.65. The reported upper bound 
on ht{mt) is valid also for the MSSM: this bound is obtained in the limit A = 0, where 
the RGE for ht is reduced to the MSSM one. Using equation (8), one can also transform 
(for fixed value of the top mass) the restriction on ht{mt) to the lower bound on tan/3. 
Clearly, this lower bound is the same both for the MSSM and for the NMSSM. 
The restrictions on A, as functions of tan/3, are presented in Fig, 1. As one can see from 
this figure, |A(mi)| < 0.7 - this bound is obtained for tan/3 > 3 and nijnt) = 0. Thus, 
the absolute upper boundary on A, derived here, is about 20% lower than IR-quasi fixed 
point A = 0.87. Such a difference arises, when the constraints on ht{mt), coming from the 
recent experimental results for the top quark mass, are taken into account. 
The behaviour of Xmax with tan f3 is nontrivial. As one can see from Fig, 1, upper bound on 
A remains almost unchanged, when 3 < tan /3 < 6, it decreases slowly, when 2 < tan /5 < 3 
and, finally, Xmax is dropped to zero for the lowest allowed values^ of tan /?. 

^In other words, Xmax is highly sensitive to ht{mt), when the latter is close to its upper boundary. On 
the other hand, the RGE's analysis approach is not able to make predictions for the couplings with such 
an accuracy. This makes the predictions for A unreliable, when the lowest values of tan/3 are considered. 
Notice however that our main results are obtained for tan/3 « 2, where such a problem does not exist. 
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Figure 1: Upper bound on |A(mt)| (Xmax) as a function of tan/? a) for K{mt) — and 
mt = leOGeV (line 1), rut = 165GeV (solid line) and rut = 170GeV (line 3), b) for 
mt = leOGeV and \K{mt) \ = 0; 0.3; 0.4; 0.5; 0.6 (lines 1,2,3,4,5 respectively) 
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It is interesting to investigate also, how the restrictions on A are affected by the experimen- 
tal error in the top mass. For this purpose I find Xmax for nit = 160GeV, = 165GeV 
and nit = ITOGeV. As one can see from Fig, la), either Ama2:(tan/3) or tan P{Xmax) vary 
only weakly, when the top mass varies form 160GeV to 170GeV. Notice however that this 
weak variation of Xmax (or tan (3) becomes nonnegligible, when being combined with the 
sensitivity of the lightest Higgs mass to the top mass via the radiative corrections [H^, 



It is clear that the difference between the pole and the on-shell top masses is nonnegli- 
gible too. This is in spite of this difference affecting only weakly on Xmax'- one obtains 
|A(mt)| < 0.65, when taking nit = 174GeV (or ht{nit) > 1/ sin/3). 

The restrictions on A, discussed above, were derived in the limit k = 0. Notice however 
that they may be considered approximately valid for \K{nit)\ < 0.3. Indeed, as one can see 
from Fig, lb), the difference between the results, derived for n{nit) = and |/t(mi)| = 0.3 
is negligible for tan/? < 2 and only weak for tan/3 > 2. More precisely, for \K{nit)\ = 0.3 
one derives |A(mt)| < 0.65 instead of |A(mi)| < 0.7. For fixed values of the top mass 
and tan j3 such a small difference does not affect the predictions for the lightest Higgs 
mass. This allows one to avoid the problems, connected with the correlations between 
the predictions for A and k. 



6. Let us return to the Higgs bosons masses problem in the MSSM and its simplest 
extension with the additional Higgs singlet. As it was discussed previously, the restric- 
tions on the lightest Higgs mass depend on the gauge coupling constants, on the Higgs 
vev's ratio tan/3, on the top quark mass, on the stops masses, which are identified here 
with the supersymmetry breaking scale, on the parameter Xt, given by eq, (3), and (in the 
NMSSM) on the coupling constant A. Here I take 1.4 < tan/3 < 6 (the lower bound on 
tan/? is derived, using the upper bound on ht{nit)), Msusy = ITeV, = and Xt = 6. 
The values of the gauge couplings and the top quark mass, as well as the restrictions on A 
were presented in previous section. The derived upper bound on the lightest Higgs mass 
ip^hmax) a function of tan/? is presented in Fig, 2. 

As one can see from this figure, the MSSM and the NMSSM upper bounds on the lightest 
Higgs mass are different in general. This difference (A/^) is maximal for tan (3 being close 
to 2. For lower values of tan/3 the difference between the models predictions decreases 
fastly to zero (together with Xmax)- For larger values of tan/? A/^ decreases too: as one 
can see from Fig, 2, the predictions of the models tend slowly to coincide for tan/? ^ 1. 
The results, presented in Fig, 2, are obtained for nit = 160GeV, = 165GeV and 
nit = 170GeV, for the cases of no-mixing (Fig, 2a) and maximal mixing (Fig, 2b). The 
largest difference between the models predictions is derived for mt = 160GeV. Thus, in 
the case of no-mixing one obtains Ah ~ 40GeV. For mt = 165GeV and mt = 170GeV it 
is obtained Ah ~ 32GeV and Ah ~ 25GeV respectively. The situation is similar for the 
case of maximal mixing. Here one obtains Ah ~ 35GeV, 28GeV and 22GeV respectively 
for mt = 160GeV, 165GeV and 170GeV. 

Thus, the difference between the models predictions varies about 15GeV, when the top 
mass varies from 160GeV to 170GeV. This is caused by different behaviour of the lightest 
Higgs mass upper bound with the top mass in the MSSM and in the NMSSM. In the 
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Figure 2: Upper bound on the lightest Higgs mass in the MSSM and the NMSSM a) 
for Xt = (no- mixing) b) for Xt = 6 (maximal mixing). The results are derived for 
rrit = IGOGeV (dashed line), rrit = 165GeV (soUd line) and for rrit = 170GeV (dotted 
hne) . 
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MSSM the experimental error ±5GeV in the top mass leads to approximately the same 
error in mh^^^. On the contrary, in the NMSSM for low values of tan/3 rrih^^^ either 
remains unchanged or decreases with the increasing of the top mass. 
There is also the difference in the behaviour of the models predictions with tan f3. As one 
can see from Fig, 2, in the MSSM rrih^^^ increases monotonously with tan/?, so that the 
absolute upper bound on is reached for tan/5 ^ 1. On the contrary, in the NMSSM 
the absolute upper bound is saturated for tan /5 = 2 — 2.5. Such a difference in the mod- 
els predictions is connected, in particular, with the difference between the pole and the 
on-shell top masses, as well as with two-loop order radiative corrections to the lightest 
Higgs mass. On the absence of these two effects the NMSSM absolute upper bound on 
rrih is saturated for tan/? ^ 1 too 0. 

It is interesting to compare the results, derived here in the case of the NMSSM, to those, 
reported in literature previously. Let us start with the case of no-mixing - in this case 
it is obtained that m/j < 120GeV. The study of this case, using the RGB's for the Higgs 
potential effective coupling constants, has been carried out in ref. - there the result 
rrih < 145GeV had been obtained. As one can see, the difference between the results is 
25GeV. This difference is caused both by stronger restrictions on A and by next-to-leading 
order corrections to the lightest Higgs mass. These two effects modify the predictions for 
rrih^^^ by 15GeV and lOGeV respectively. 

In the case of maximal mixing it is obtained rrih < (130 — 135) GeV (depending on the 
top mass). This result is compared to those, reported in ref's [|l^, |l^, ^ - there the 
upper bound ruh < 155GeV is obtained. Again, the difference between the results is 
(20-25) GeV. However, unlike the case of no- mixing, the dominant contribution to this 
difference (15GeV) comes now from two-loop order radiative corrections. 

7. Thus, the problem of the restrictions on the lightest Higgs boson mass in the next-to- 
minimal supersymmetric standard model has been analyzed. To obtain these restrictions, 
the constraints on coupling constants of theory were taken into account. These constraints 
come from the renormalization group equations analysis and from the recent experimental 
results for the top mass. The experimental constraints on the top mass are very important: 
they make the restrictions on coupling constants much stronger. The another important 
effect is connected with two-loop order radiative corrections to the lightest Higgs mass. 
Due to these two effects the derived NMSSM upper bound on the lightest Higgs mass is 
about (20-25) GeV lower than the one, reported in literature previously. 
There is a large difference between the MSSM and the NMSSM constraints on the lightest 
Higgs mass. Namely, the upper bound on the lightest Higgs mass in the NMSSM can 
be about 40GeV larger than in the MSSM. Such a large difference occurs for tan/3 being 
close to 2. Due to this difference the NMSSM with the low tan/? scenario is, generally 
speaking, far out from being excluded. Notice however that the NMSSM upper bound on 
the lightest Higgs mass is saturated in the narrow region of parameters space, where the 
singlet Higgs sector is decoupled and the detectable Higgs bosons masses are described 
by the MSSM-type mass matrices. To compare the models predictions in wider region of 
parameters space, the predictions for the supersymmetry breaking parameters of theories 
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are necessary. 
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APPENDIX 

At tree level the 3x3 symmetric mass matrix for the fields $i, $2 and Ni is the following: 

M|^^ = M| cos^ 2/? + A V sin^ 2/3 
M| = -(M| - AV)sin2/?cos2/3 



= 2X^37] + {Ax + 2KV3)Xvsm2f3 

sin 2(3 

r2 



■513 

Ml, = iMl-XY)s\n'2(3-2Xvs ''t!n7 (A-l) 



Ms,, = -iAx + 2KVs)X7]cos2l3 
-2 

■533 



Ml, = ^K^vl + KV3A^- 



AxXt]"^ sin 2(3 

2^3 



where Ax, A^, are trilinear SUSY breaking parameters of the Higgs potential |T0|, |TT|, [13 
IH], is the singlet vev. The fields $1, $2, ^1 are given by equation (4). 
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